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Abstract. One of the goals of the ground-based support program for the COROT and MONS/R0mer satellite 
missions is to select and characterise suitable target stars for the part of the missions dedicated to asteroseismology. 
While the global atmospheric parameters may be determined with good accuracy from the Stromgren indices, 
careful abundance analysis must be made for the proposed main targets. This is a time consuming process 
considering the long list of primary and secondary targets. We have therefore developed new software called VWA 
for this task. The VWA automatically selects the least blended lines from the atomic line database vald, and 
consequently adjusts the abundance in order to find the best match between the calculated and observed spectra. 
The variability of HD 49434 was discovered as part of COROT ground-based support observations. Here we 
present a detailed abundance analysis of HD 49434 using VWA. For most elements we find abundances somewhat 
below the Solar values, in particular we find [Fe/H] = —0.13 ± 0.14. We also present the results from the study 
of the variability that is seen in spectroscopic and photometric time series observations. From the characteristics 
of the variation seen in photometry and in the line profiles we propose that HD 49434 is a variable star of the 7 
Doradus type. 

Key words, methods: data analysis - techniques: spectroscopic stars: abundances - stars: fundamental parameters 
- stars: individual: HD 49434 - stars: variables: 7 Dor 



1. Introduction 

We first observed HD 49434 (F2V-type) with the ELODIE 
spectrograph in December 1998 as part of a program with 
the aim of characterising and selecting potential targets 
for the asteroseismology space mission COROT (Catala 
2001). The strategy is to obtain a number of high resolu- 
tion, high signal-to-noise spectra of each potential target. 
For any stars that present line asymmetries or variable 
line profiles we then carry out more complete monitoring. 

Before these observations the star was not consid- 
ered to be variable, but we discovered systematic low- 
Send offprint requests to: H. Bruntt, e-mail: 
bruntt(5ifa.au.dk 

* Based on observations obtained at Observatoire d'Haute 
Provence, France and at the Observatory of Sierra Nevada, 
Granada, Spain 



amplitude variations of the line profiles. Since then we 
have carried out more extensive spectroscopic and photo- 
metric observations. 

The purpose of the present paper is to derive the global 
atmospheric parameters of HD 49434 and to perform an 
analysis of its photometric and spectroscopic variations. 

The French COROT as well as the Danish 
MONS/R0mer satellite are both asteroseismology mis- 
sions which will study mostly solar-like stars. Before the 
launch of these missions (around 2004-5) we need an accu- 
rate estimation of the atmospheric parameters and abun- 
dances of the proposed target stars. First of all stars with 
certain parameters may not be suitable for a detailed as- 
teroseismological study, eg. chemically peculiar stars or 
stars with high rotation rate. Also, in order to carry out 
the detailed asteroseismological modelling of the stars we 
need to be able to constrain as many of the global param- 
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eters of the stars as possible. Around 10 and 25 primary 
targets will be chosen for COROT and MONS/R0mer, re- 
spectively. It is a very time-consuming process to carry out 
the abundance analysis of all potential targets manually 
and we have therefore developed a fast semi-automatic 
abundance analysis procedure called vwa, which will be 
described here. The result for HD 49434 is presented, and 
we will also make comparisons with stars for which abun- 
dance analyses have already been carried out manually. 

In Section 2 we describe the spectroscopic observa- 
tions and in Section 3 we describe the VWA software and 
the result of the abundance analysis. In Section 4 we de- 
scribe the analysis of the photometric time series, while 
in Section 5 we analyse the line profile variations in the 
spectroscopic time series. 



2. Spectroscopic observations 

HD 49434 is an F2V-type star of magnitude V = 5.74 
(cf. Table 2). We monitored the star spectroscopically 
from Observatoire de Haute-Provencc, using the ELODIE 
spectrograph, attached to the 1.93 m telescope. The spec- 
trograph is a cross-dispersed spectrograph, providing a 
complete spectral coverage of the 3 800-6800 A region, 
with a resolving power of about R = 45 000 (Baranne et 
al. 1996). 

The first spectroscopic observations of HD 49434 were 
carried out in December 1998 as part of the search for po- 
tential COROT targets. The spectra revealed the presence 
of variable asymmetries of low amplitude in the line pro- 
files, indicative of non-radial pulsations. Hence, the star 
was monitored more thoroughly in subsequent observing 
runs. 

On January 26, 2000, HD 49434 was monitored con- 
tinuously for 90 minutes, with one 5 min spectrum every 
7 minutes. The weather and seeing conditions were very 
good, and the average S/N ratio at 6 000 A is 250 per 
velocity bin of 10 kms -1 . 

The star was re-observed on December 10, 2000, al- 
most continuously for 320 minutes, during which as many 
as 51 five-minute spectra were recorded. The observing 
conditions of this second monitoring were not as good as 
the previous one, so that the resulting average S/N ratio 
at 6 000 A is only about 170 per pixel per velocity bin of 
10 kms" 1 . 

The data used for the abundance analysis of HD 49434 
are the spectra obtained during the best weather condi- 
tions: 12 spectra from January 26, 2000. On the other 
hand, all of the observed spectra were used for the study 
of the variability of line profiles. Since the data analysis of 
the spectra was made by two groups independently, we will 
describe the data reduction in separate sections, i.e. abun- 
dance analysis in Section 3.1 below and the variability of 
line profiles in Section 5.1. 



3. Abundance analysis 

In this section we will describe the abundance analysis 
of HD 49434. For this purpose we have developed new 
software for the semi-automatic analysis of spectra. The 
abbreviation for the software is vwa which explains what 
the program can determine from the spectrum, i.e. vsuii 
of the star, wavelength shift (due to radial velocity), and 
abundance analysis. 

3.1. Data reduction 

The bias subtraction, flat fielding, scattered light subtrac- 
tion, and extraction of the spectral orders were done with 
iraf. The wavelength calibration was made using a com- 
bined Th-Ar comparison spectrum. We made low-order 
cubic spline fits to make the best continuum estimate. We 
made sure that the overlap between orders was better than 
0.5%. We note that we have not used the region near the 
hydrogen lines for abundance analysis. The hydrogen lines 
cover several spectral orders and thus a determination of 
the continuum level is difficult. One possibility to over- 
come this problem is to use the average of the continuum 
level for orders adjacent to the hydrogen lines. This has 
been done by Lastennet et al. (2001) and they have esti- 
mated the effective temperature from the Ha line to be 
T cff = 7250 ± 250 K. 

The signal to noise of the combined spectrum (12 spec- 
tra) at the maximum of a spectral order around 6 000 A is 
S/N = 400 per velocity bin of 10 km s" 1 . 

3.2. Automatic abundance analysis with vwa 

The vwa package can perform three different tasks: 
Automatic selection of lines for the determination of 
usini, determination of any systematic shift of lines (due 
to radial velocity), and abundance analysis the last being 
the most extensive part of the program. 

The first two tasks are not fully implemented in vwa, 
but v sin i and v ra( j can be determined quite easily by com- 
parison of the observed spectrum and several synthetic 
lines: We convolve the synthetic spectrum with the in- 
strumental profile and different rotation profiles and deter- 
mine the best fit (we use the program rotate by Piskunov 
1992). We find vs'mi = 84 ± 4 kms -1 which is in agree- 
ment with Lastennet et al. (2001) who found vsini = 79 
km s -1 from the same data set. From the wavelength shifts 
of non-blended lines we determine the helio-centric radial 
velocity to be w ra d = —12.3 ± 1.0 kms -1 . Gernier et al. 
(1999) found w ra d = —14.0 ± 0.8 kms -1 from Hipparcos 
data. We remove the wavelength shift caused by the ra- 
dial velocity, thus the wavelength scale is the same as the 
laboratory system. 

3.2.1. Input: model atmosphere and atomic data 

To use VWA one must calculate an appropriate model for 
the atmosphere of the star. We have used a modified ver- 
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Fig. 1. Observed (thin line) and synthetic (thick line) spectrum around twelve lines which were used for the abundance deter- 
mination. The lower portion of the each plotted line shows the relative difference between the observed and synthetic spectrum. 
The name of the element and the central wavelength (in Angstr0m) of each fitted line is given at the top of each plot 
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sion of the ATLAS9 code (Kurucz 1993) as described by 
Kupka (1996) - see also Smalley & Kupka (1997). In this 
version of ATLAS9 the turbulent convection theory in the 
formulation of Canuto & Mazzitclli (1991, 1992) is imple- 
mented. 

To calculate the model one needs to know the basic 
parameters of the star, i.e. T e s,\ogg, and [Fe/H]. These 
parameters may be estimated from the Stromgren indices 
of the star by using the program TEMPLOGG (Rogers 1995) 
which automatically chooses the most appropriate among 
several published photometric calibrations. To determine 
the basic atmospheric parameters of the star TEMPLOGG 
then interpolates the de-reddened photometric indices in 
the atmospheric model grid by Kurucz (1979) with the 
improvements suggested by Napiwotzki et al. (1993). 

For the abundance determination we make use of 
a complete list of atomic line data for the entire ob- 
served spectral range. The line list is extracted from the 
VALD database (Kupka et al. 1999, Ryabchikova et al. 
1999) where all lines deeper than 1% of the continuum 
are included. From the line list (about 13 000 lines for 
HD 49434) the least blended lines will be selected by vwa. 



The selection consists of two steps: 1) Selection based 
on the analysis of the degree of blending based on the 
information in the line list from the VALD database. 2) 
These lines are then analysed in the observed spectrum 
and any "suspicious" lines are rejected. These steps are 
described in the following two sections. 

3.2.2. Selecting lines: atomic line list 

For each atomic line extracted from the VALD database we 
examine the line depths of the neighbouring lines. For a 
given line we examine lines within a certain range specified 
by the user; for HD 49434 we have used 3.0xFWHM of an 
unblended line (= 3.0 x 2.2 A= 6.6 A). The line depths of 
the neighbouring lines are convolved by a Gaussian with 
a FWHM of 1.5xFWHM of an unblended line, i.e. the 
line depths of lines far away from the central line are sup- 
pressed. The selection of lines from the VALD line list is 
based on the degree of blending which is described by 
three parameters. The user may set his own limits for these 
parameters, and thus determine which lines are selected. 
The parameters are: 1) The depth ratio of the central line 
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to the deepest (non-convolved) line. This value must be 
greater than 1 but can specified by the user (we used 1.35). 
2) The number of neighbours with a (convolved) line depth 
of 25% of the central line must be lower than some limit 
(we used 5). 3) The sum of the (convolved) line depths 
must be lower than some number (we used 0.7). 

3.2.3. Rejecting lines: the observed spectrum 

When the lines have been selected from the list of atomic 
lines vwa analyses the observed spectrum in the region 
around each selected line. The aim is to automatically 
reject lines that are not suitable for abundance analysis. 
This is done by fitting a Voigt profile to the observed spec- 
trum. A number of "selection parameters" specified by the 
user determines if a line will be used: the minimum depth 
of the fitted profile, the maximum and minimum width of 
the profile, the maximum value of x 2 for the profile fit, 
and the maximum flux asymmetry between the left and 
right part of the observed line. 

Optionally the user is presented with a plot of the ob- 
served spectrum in a region around each of the selected 
and rejected lines. The user may then interactively choose 
if a rejected (accepted) line should be included (rejected) 
in the abundance analysis. In practice this part of the VWA 
program creates two lists of the lines that will be rejected / 
accepted in the final selection of lines to be used for abun- 
dance analysis, i.e. these line lists over-rule the automatic 
selection. The line lists may also be edited manually if, for 
example, the user already knows which lines he wants to 
use for the abundance analysis (in accepted line list) - or 
if the user knows that the atomic parameters of certain 
lines are known to be uncertain (in rejected line list). 

3.2.4. Two methods for abundance determination 

The abundance determination is made for each of the se- 
lected lines. It can be done by either of two methods: 1) 
By the measurement of equivalent widths of lines in the 
observed spectrum which are then given as input to the 
width9 code by Kurucz (1993). 2) By the adjustment of 
the abundance of the element of the central line of the 
blend until the computed and observed spectrum around 
each line match, i.e. they have the same equivalent width. 
For the computation of the synthetic spectrum we used the 
program Synth (version 2.5) by N. Piskunov (see Valcnti 
& Piskunov 1996). 

The first method can only be used for stars where 
blending is not a problem, i.e. slowly rotating stars 
(wsini < 15 kms -1 ). The second method is much slower 
but it is the only option for fast rotators like HD 49434. 

3.2.5. Deriving the abundances 

For this study we have used the fitting of the abundance 
in the synthetic spectrum in a region 4 A on both sides 
of the central line (the typical FWHM of a single line is 



about 2.2 A at 6 000 A). After all lines have been fitted 
the user is presented with a plot of the region around each 
line in which both the synthetic and observed spectrum is 
shown. Also, all lines with a line depth deeper than 25% of 
the central line are clearly marked. Several fitted lines are 
shown in Figure 1 . In this "visual inspection" part of VWA 
the user can click on different sections of the plot which 
determines if the user thinks the line is fitted correctly, 
if there seems to be a problem with the continuum level, 
if the line is severely blended, if the line should simply 
be ignored, or if the automatic fit has failed. In the last 
case these lines can be fitted interactively by a separate 
program (vwa-interact) . 

The abundance of each element is then calculated us- 
ing only the lines accepted by the user. Additional con- 
straints on which lines to use may be imposed here, i.e. to 
use only the least blended lines. For this the "sensitivity" 
parameter of the line can be used to select the best lines 
or to give higher weight to certain lines. We define the 
sensitivity parameter as the change in equivalent width 
when increasing the abundance, i.e. for a non-blended line 
this would be the slope of the "curve of growth" near the 
abundance of the fitted line. The parameter is measured 
for each line during the automatic determination of the 
abundance and we find that it depends both on the prop- 
erties of the line and the degree of blending. 

To be specific, the sensitivity, Si, for each line i, is de- 
fined as S t = A(EQW)/A(log N/N tot ). Typical values lie 
in the range 0.3-1.2 mA/ dex. For the present data we as- 
sume that the measured equivalent widths of all lines are 
measured with an accuracy of 5%, hence the error estimate 
of the abundance for a given line lies in the range 0.17-0.04 
dex. To get a realistic error estimate for each line one must 
include systematic errors due to eg. erroneous log gf val- 
ues (i.e. oscillator strengths), wrong determination of the 
continuum level, the possibility of erroneous abundances 
(or log gf values) of the lines that blend with the line, and 
uncertain stellar parameters used for the computation of 
the atmosphere model. 

When calculating the abundance of each clement based 
on several lines we assign weights to each line which de- 
pend on the sensitivity of the line. When computing the 
weights we have added an additional error of 0.07 dex 
which is the contribution from systematic errors as was 
discussed above. For the final abundance result of each ele- 
ment we calculate the weighted mean, where the weights 
are given by Wf 1 = (0.05/S* 4 ) 2 + (0.07) 2 . The standard 
error on the derived abundance is also calculated using 
these weights. 

3.2.6. Reliability of vwa 

We have tested the program on two stars for which careful 
analyses have already been made: The Sun (Grevesse & 
Sauval 1998) and FG Vir (A-typc star, usini = 21.3 ± 1.0 
kms -1 ; see Mittermayer 2001). For these stars we find 
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Table 1. Abundances of elements in HD 49434. The weighted mean abundance is given as the logarithm of number of atoms 
of a given element to the total number of atoms. The weighted standard deviation of the mean is given in parenthesis, n is the 
number of lines that were used, while the last two columns are the abundances measured for the Sun (Grevesse & Sauval 1998) 
and the difference between the Sun and HD 49434 



Element 


Ion 


logAVAU 


n 


(logiV/iV tot ) G 


A 


C 


I 


-3.46(0.13) 


1 


-3.52 


+0.06 


Mg 


I 


-4.55(0.07) 


3 


-4.46 


-0.09 


Si 


I 


-4.64(0.08) 


2 


-4.49 


-0.15 


Si 


II 


-4.57(0.06) 


3 


-4.49 


-0.08 


Ca 


I 


-5.51(0.06) 


7 


-5.68 


+0.17 


Sc 


II 


-8.89(0.11) 


2 


-8.87 


-0.02 


Ti 


II 


-7.16(0.07) 


2 


-7.02 


-0.14 


Cr 


II 


-6.45(0.05) 


4 


-6.37 


-0.08 


Fe 


I 


-4.67(0.02) 


16 


-4.54 


-0.13 


Fe 


II 


-4.43(0.07) 


2 


-4.54 


+0.11 


Ba 


II 


-9.75(0.09) 


3 


-9.91 


+0.16 



good agreement with the previously derived abundances, 
i.e. within 0.1 dex for all elements. 

We note that two future papers which make use of 
VWA are in preparation, i.e. Rasmussen et al. (2002) and 
Bikmacv et al. (2002). In the later paper we study two 
A-type stars and both have low vsin? (9 and 21 kins" 1 ). 
Thus we may compare the result of VWA with a "classi- 
cal" abundance determination, i.e. by measuring equiva- 
lent widths of non-blended lines. The two methods agree 
within 0.03 dex for all elements. We conclude that VWA is 
indeed a reliable tool for fast semi-automatic abundance 
analysis. 

3.3. Abundances for HD 49434 

We have estimated the fundamental atmospheric parame- 
ters of HD 49434 from the Stromgrcn photometry given in 
Table 2. We have used the templogg code (Rogers 1995, 
see also Kupka & Bruntt 2001) and we find the atmo- 
spheric parameters of HD 49434 to be T off = 7 300+200 K, 
\ogg = 4.14 + 0.20, and [Fe/H]= -0.01 + 0.20. The quoted 
errors are approximate and are dominated by the uncer- 
tainty from the photometric calibrations. 

We determined the microturbulcncc parameter to be 
£t = 2.1 ± 0.5 kms -1 (cf. discussion in Section 3.3.1). We 
used macroturbulcncc f maC ro — 5.0 kms -1 and vsini = 
84 + 4 kms" 1 . The choice of w maC ro is not important here, 
since the width of the lines is dominated by the rotational 
broadening. 

Several lines fitted by VWA are shown in Figure 1. We 
have also shown the relative difference between the syn- 
thetic and observed spectrum. Neighbouring lines with a 
line depth of at least 25% of the fitted line are also shown. 
The length of the tick marks show the relative line depth 
of each neighbouring line compared to the fitted line. 

From Figure 1 it can be seen that several lines that 
were used for the abundance analysis are blended. Due to 
the high v sin i of HD 49434 we simply had to use the least 
blended lines. Our strategy was to determine the Fe abun- 
dance accurately from non-blended lines. Consequently, 



Table 2. Stromgren photometric indices taken from Hauck & 
Mermilliod (1998) for the four stars used for the HD 49434 
observations from Sierra Nevada 



Star 


V 


b-y 


mi 


Cl 


H/3 


HD 48922 


6.77 


-0.015 


0.142 


0.940 




HD 49933 


5.78 


0.270 


0.127 


0.460 


2.662 


HD 50747 


5.45 


0.095 


0.172 


1.131 


2.833 


HD 49434 


5.74 


0.178 


0.178 


0.717 


2.755 



we held the abundance of Fe fixed, and allowed VWA to 
select lines of other elements that were mildly blended by 
Fe lines, eg. the Ca lines at 5594.5 and 6102.7 A. More 
extreme examples are the C, Ni, and Sc lines at 5052.2, 
5614.8, and 5526.8 A where the lines are heavily blended 
and the abundance determination is unreliable, i.e. the Si 
are 0.49, 0.24, and 0.22. If we again assume that the equiv- 
alent width is known to 5% the corresponding error is 0.1 
dex for the C lines and 0.2 dex for the Ni and Sc lines. 
In addition, it is difficult to establish the location of the 
continuum level because the lines are so broad. Realistic 
error estimates lie in the range 0.2-0.4 dex for these three 
lines. Note that we did not find enough lines to give a re- 
liable estimate for the abundance Ni and for C we found 
only one line that was usable. 

The result of the abundance analysis is given in Table 1 
and for comparison we also give the solar abundances from 
Grevesse & Sauval (1998). The abundances and errors are 
calculated using weights that are estimated from the er- 
rors on the individual lines (cf. Section 3.2.5). The error 
estimates that we give in Tabic 1 do not include the sys- 
tematic error due to the uncertainty of the atmospheric 
model parameters which will be discussed in Section 3.3.1. 

Generally we find that the abundances of most ele- 
ments lie around the Solar value within the error bars. 
Only for Fe and Ca do we find enough lines (16 and 7) to 
make a truly reliable estimate of the abundances. 

The Fe abundance is log N Fc /N tot = -4.67 + 0.02 from 
16 Fel lines (we quote the internal weighted error; the 
internal RMS scatter is 0.18 dex). If we assume that the 
photospheric hydrogen and helium abundance is the same 
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Fig. 2. Differential magnitudes of HD 49434 minus comparison star HD 48922 in the Stromgren b band (asterisks). The small 
points are magnitude differences between the two comparison stars HD 49933 and HD 48922 (shifted to allow comparison). 
Individual days are marked in every panel as the corresponding HJD 
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in HD 49434 and the Sun this abundance corresponds to 
[Fe/H] = -0.13 ±0.14 (including systematic errors). This 
agrees roughly with the result from the Stromgren mi 
index which yields [Fe/H] = -0.01 ± 0.20. 

An interesting result is the ratio [Ca/Fe] = +0.30 ± 
0.21 which is a quite high value for a star with Solar-like 
abundances. 

We note that we have adjusted the log gf for the Can 
line at 5857.451 A line by comparing with the solar spec- 
trum, i.e. from log gf = 0.257 to 0.550. The log<?/ val- 
ues for the Si lines are determined solely from theoretical 
calculations. Therefore we have computed the synthetic 
spectrum for the Si-lines we have used - using the solar 
abundance - and comparing this with the observed solar 
spectrum. We find that the log gf values seem to be right 
for the lines we have used. 

Lastennet et al. (2001) have determined fundamental 
atmospheric parameters of HD 49434 by comparison of 
the observed and synthetic photometric colours. The star 
is among the nine potential COROT targets stars they 
have studied. From the combination (B — V), (U — B), 
(6 - y) they find log g = 4.0 ± 0.4, T cff = 7 240 ± 100 K 
and [Fc/H]= -0.1 ± 0.2 (from Table 2 and Figure 7 and 8 
in Lastennet et al. 2001). Thus, the metallicity they find 
is in agreement with our result from spectroscopy and all 
the fundamental parameters are in agreement with the 
estimates based on Stromgren photometry 

We finally note that HD 49434 is included in the Ao 
photometry catalogue of bright B and A stars by Vogt 
et al. (1998) but they find no evidence for the star being 
chemically peculiar. 



3.3.1. Accuracy of the derived abundances 

The errors quoted in Table 1 for the abundances are the 
weighted internal errors (cf. Section 3.3.1). Here we will 
discuss the error contribution from the uncertainty of the 
atmospheric model parameters. The atmospheric model 
we used for the final abundance result have the param- 
eters T e g = 7 300 K and log g = 4.1. To investigate the 
effect of the choice of atmospheric model we calculated 
Fe abundances for models with T cS = 7 100, 7 300, 7 500 K 
and \ogg = 4.1 and one model with T e g = 7 300 and 
log g = 4.3. We find that the effect of increasing \ogg by 

0. 2 dex decrease the Fe abundance by about 0.04±0.05 dex 
(sd. of mean) which is not significant. On the other hand, 
if we increase the temperature of the atmospheric model 
by 200 K we find an increase of Fe by 0.10 ± 0.05 dex, 

1. e. a significant effect. 

Another input parameter for the calculation of the syn- 
thetic spectrum is the microturbulence. The value used in 
the calculation of the Kurucz atmosphere model does not 
affect the derived abundances significantly, but when cal- 
culating the synthetic spectrum the effect is indeed signif- 
icant. When changing the microturbulence the change in 
equivalent width depends on the strength of the line. By 
using this fact, one can adjust the microturbulence until 
the abundance of individual lines do not correlate with 
equivalent width. We could not do this very accurately 
due to the high rotational velocity of HD 49434. Thus, we 
estimate the microturbulence to be £ t — 2.1 ± 0.5 kms -1 . 

We find that when increasing £ t by 0.5 kms -1 the 
abundance of Fe decreases by 0.10 ± 0.04 dex, i.e. a signif- 
icant contribution to the uncertainty of the abundance. 
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Fig. 3. Power spectra corresponding to the variable star 
HD 49434 minus the comparison star HD 48922, in the up- 
per panel, and that corresponding to the comparison stars 
HD 49933 minus HD 48922 in the bottom panel. The scales 
are milli-magnitudes and cycles/day 
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Fig. 4. Location of HD 49434 (star symbol) in the HR diagram 
together with the known 8 Scuti (circles, Rodriguez & Breger 
2001) and bona fide 7 Dor (triangles, Handler & Kaye 2001) 
pulsators. The edges of the 7 Dor region around (b — y)o = 0.2 
are from Handler (1999). The absolute magnitude of HD 49434 
is My = 2.73 ±0.07 and is derived from its Hipparcos parallax 
(ESA 1997) 
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We conclude that the contribution to the error of abun- 
dances due to uncertain model atmosphere parameters, 
i.e. T c ff and microturbulence, is of the order 0.14 dex. 
For a realistic estimate of the abundance of elements in 
HD 49434 this contribution must be added to the internal 
errors (sd. of mean) given in Table 1. For example for the 
metallicity of HD 49434 we get [Fe/H] = -0.13 ± 0.14, 
i.e. the uncertainty of the model parameters is the main 
error source. 

4. Stromgren photometry time series 

We have observed HD 49434 with a Stromgren photometer 
on the 0.9m telescope at Sierra Nevada near Granada, 
Spain. We used a 45" diaphragm to observe HD 49434 
and the comparison stars HD 48922=C1, HD 49933=C2, 
and HD 50747=03. The star HD 49933 is a V = 5.78 
double system with a faint companion of magnitude V 
= 11. This last star was included within the diaphragm 
because of its proximity but had no influence on the final 
differential photometry. The Stromgren parameters for the 
observed stars are given in Table 2. 

Continuous photometric time series of comparison 
stars HD 48922, HD 49933 and HD 49434 were made in 
order to remove the sky transparency fluctuations to be 
able to see the variations of the target star. Magnitude 
differences between HD 49933 and HD 48922 are constant 
within 3.7 mmag (rms) and clear intrinsic variations can 
be seen for HD 49434. The light curve is shown in Figure 2. 

The time series consist of data from nine nights 
spread unevenly over one month starting on July 2001. 
Unfortunately the weather conditions were not good dur- 
ing our first night of observation and HD 49434 was ob- 
served again a week later. Observations were possible only 
three weeks later where the star was included in a program 
where other targets were also measured. Thus the dedi- 



cated observing time was limited to only a few measure- 
ments but the data could be used to confirm the long-term 
variations that were seen in the first days. 

In total the time series consists of 295 data points in 
the four Stromgren bands (uvby) . In Figure 2 we show the 
time series of HD 49434 - HD 48922 (target - comparison 
star 1) and HD 49933 - HD 48922 (comparison star 2 - 
comparison star 1). The magnitude differences are plotted 
for the Stromgren b band. 

From Figure 2 it is clear that the intrinsic photometric 
variations are associated with HD 49434. At first sight the 
star does not seem to have periodic variations. A power 
spectrum was calculated in order to search for possible 
frequency components which is shown in Figure 3. It can 
be seen that with the present data, no clear discrete fre- 
quency components are present. Instead the low frequency 
region has a significant contribution for HD 49434 which 
is not seen for the comparison stars. 

This behaviour is typical among the recently discov- 
ered 7 Dor stars (Zerbi et al. 1999) which are g-mode 
pulsators. We believe that HD 49434 is a variable star of 
this type: the star is classified as a main sequence F2- 
type star from Stromgren photometry and its position in 
the Hertzsprung-Russell diagram in Figure 4 confirm our 
hypothesis. Because we have not detected any clear pe- 
riodicity it is not possible to test the phase differences 
among the different photometric colours which is a clear 
signature of the non-radial pulsating nature of the 7 Dor 
stars - see Garrido (2000) for details. 

5. Spectroscopic monitoring of HD 49434 

In the following we will discuss the analysis of the time 
series of spectra of HD 49434 to look for any systematic 
variation of line profiles. 
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Fig. 5. The residual of the mean profile for the January 2000 series. Left: in velocity - time space; the maximum amplitude 
is of the order of 3 x 10~ 3 ; Right: Spectral power density of the same data, in apparent order - frequency space, after the 2D 
Fourier transform 
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Fig. 6. The same as Figure 5 but for the December 2000 series; the maximum amplitude of the residuals is of the order of 
4 x 10~ 3 



5.1. Data reduction 

All data were reduced on-site, using the automatic 
INTER-TACOS procedure (Baranne et al. 1996), result- 
ing in wavelength-calibrated, flat-field corrected spectra. 
Each spectrum is segmented into 67 segments correspond- 
ing to the 67 orders of the spectrograph present on the 
detector. We then performed a large scale normalisation 
of the overall spectrum, so that the resulting spectra are 
normalised to unity all across the covered wavelength do- 
main, even though some local problems still remain with 
this normalisation on the edges of the orders. 



The ELODIE spectrograph provides a wide enough 
spectral coverage to give access to a large number of 
photospheric lines. In order to take advantage of this, 
and to fully exploit the information contained in these 
numerous photospheric lines, we used the Least-Square 
Deconvolution (LSD) technique (Donati et al. 1997) to 
analyse the variations of a "mean" photospheric line. In 
this method, a line pattern function is constructed, con- 
taining all the lines supposedly present in the spectrum 
as Dirac functions, with heights set to the central line 
depths as calculated by Kurucz's (1979) "SYNTHE" pro- 
gram. The observed spectrum is then deconvolved with 
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this line pattern function, yielding a "mean" photo-spheric 
line profile. With this technique, line blends are automat- 
ically taken into account when all the lines present in the 
spectrum are considered. Note that the depth of the result- 
ing "mean" line has no physical meaning, but that time 
variations of this depth, as well as line profiles can be ac- 
curately analysed with this technique. We used a Kurucz 
model for T c g = 7 000 K and log 5 = 4, sufficiently close 
to the values derived in this paper, for constructing the 
line pattern function. 

A total of 1700 lines were finally used in this analysis. 
In the following of this work, we shall call "mean" photo- 
spheric line, the line constructed with the LSD technique, 
while the adjective "average" will be reserved for time av- 
erages of line profiles. 

All mean profiles in each series were normalised to the 
same equivalent width. This procedure eliminates the ef- 
fects of the large scale variations of the mean line pro- 
file, due for instance to normalisation problems that may 
vary from spectrum to spectrum. The time average of the 
mean profiles was then constructed for each series, and 
subtracted from the mean line profiles of the series. 

We then performed a two-dimensional Fourier trans- 
form analysis, following the method described in Kcnnclly 
et al. (1996). This method is well suited for the analysis of 
the spectroscopic signature of non-radial pulsations. Each 
line profile was interpolated on a grid representing stellar 
longitudes, transforming velocities across the line profile 
into longitudes on the stellar equator using the relation 
v — v r = vsini sin 0, where v is the velocity in the pro- 
file, v r is the radial velocity of the star, measured to be 
Vrad — —12.3 ± 1.0 kms -1 , <f) is the corresponding stellar 
longitude, and vsini is the projected rotation velocity of 
the star, estimated to be 84 ± 4 kms^ 1 . 

For the time series from January 2000, the spectra were 
recorded at constant time intervals during the whole moni- 
toring, while a few short gaps exist in the December 2000 
series. For this latter series, prior to the 2D Fourier trans- 
form, we had to build up a grid of constant time intervals, 
filling the gaps during which no spectrum was obtained 
with null data. 

After these steps, the data were submitted to a 2D 
Fourier transform, resulting in a representation of the data 
in the (frequency - apparent m) space. The apparent m is 
related to the structure of the modes present at the stellar 
surface, without being identical to the usual azimuthal 
order m, since the apparent m depends for instance on 
the inclination angle i of the rotation axis of the star with 
respect to the line of sight. The apparent m is nevertheless 
indicative of the surface structure of the modes detected. 

5.2. Detected variability of spectral lines 

Figure 5 (and Figure 6) presents the results of this anal- 
ysis for the scries from January 2000 (and December 
2000). At both epochs we detect the presence of line pro- 
file perturbations, crossing the profile from blue to red. 



Such perturbations arc most probably the signature of 
non-radial pulsation modes at the stellar surface. The 
Figures 5 and 6 suggest the presence of one or several 
high-degree modes. These modes are particularly obvious 
in the January 2000 series, owing to the high signal-to- 
noise ratio of the spectra, while they are less conspicu- 
ous, although present without doubt, in the lower quality 
December 2000 data. The amplitudes of the residuals are 
of the order of 2-3 x 10 -3 at both epochs. They can be eas- 
ily detected thanks to the very high signal-to-noise ratio 
(S/N) of the LSD mean line profiles combining hundreds 
of spectral lines, which reach S/N=3 000 in the December 
2000 series and S/N=6 000 in the January 2000 series. 

These modes seem to be long-period modes. This is 
particularly clear in Figure 6 (right), where the frequency 
of the detected modes are at the resolution limit of our 
dataset, i.e. 0.1 mHz, indicating a period longer than 160 
minutes, that is to say longer than half the time span of 
the monitoring. These results therefore confirm those of 
the photometric observations, and suggest that HD 49434 
is probably a variable of the 7 Doradus type. 

The modes we detect in the January 2000 series are 
high-degree modes, with an apparent order of 6. In the 
December 2000 series, our lower quality data still indicates 
non-radial modes, but with more power on the apparent 
order 2, some signal still appearing at apparent order 8. 

6. Conclusions and future prospects 

We have presented the detailed abundance analysis of 
HD 49434. The analysis could not be made by sim- 
ple equivalent width measurements since most lines are 
blended, because the star has high projected rotational 
velocity. We have described the vwa software which se- 
lects the least blended lines and derive abundances semi- 
automatically. 

We find that the mctallicity of HD 49434 is some- 
what below the Solar value, i.e. [Fe/H]= -0.13 ± 0.14. 
The quoted error includes the estimated error due to the 
uncertainty of the atmosphere parameters, especially T e g 
and microturbulence. 

The metallicity we find agrees with the results from 
Stromgren photometry and the work by Lastennet et al. 
(2001) who have compared observed and synthetic photo- 
metric colours. When taking into account the uncertainty 
of the basic atmosphere parameters the accuracy of VWA 
is only slightly more accurate than the photometric meth- 
ods. The advantage of analysing the spectra with VWA is 
that we obtain estimates for several elements. 

We have already used vwa for the analysis of several 
other stars which will be published soon (Bikmaev et al. 
2002, Rasmussen et al. 2002). In the paper by Bikmaev 
et al. we have analysed two A-type stars, which have low 
vsini. Thus we have made a detailed comparison of the 
result obtained with VWA and with the "classical" method 
based on equivalent width measurements of non-blended 
lines the agreement is excellent for nearly all lines (abun- 
dances < 0.03 dex). 
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With the confirmation that vwa can be used for reli- 
able semi-automatic abundance analysis we will proceed 
to use the software for analysis of the proposed main tar- 
gets for the asteroseismology satellite missions COROT 
and MONS/R0mer: the spectra have already been ob- 
tained. 

For our future use of vwa we expect that for stars with 
low vaini (< 40-50 kins" 1 ) it will be possible to con- 
strain the atmospheric parameters (eg. microturbulence 
and T e g) by adjusting these so the abundances found from 
lines of the same element do not correlate with the line pa- 
rameters (eg. equivalent width and excitation potential). 
For stars with higher i>sini (> 100 kms -1 ) the selection 
of atomic lines for abundance analysis becomes exceed- 
ingly difficult. For such stars the metallicity and the basic 
parameters must be determined mainly from photometric 
colours, but it is still important to obtain spectra to be 
able to put constraints on T e g (from the hydrogen lines) 
and v sin i. For example, stars with too high v sin i may not 
be suitable for seismology since this will make the analysis 
much more complex. 

In the present paper we have also presented the photo- 
metric light curve of HD 49434 based on observations dur- 
ing part of nine nights over a period of one month. We find 
that the star is indeed variable with excess power around 
1-5 cycles per day, but no periodic signal could be found 
from the limited data set. We have also analysed the spec- 
troscopic line profiles which show evidence for long-period 
variability at low amplitude. This behaviour is characteris- 
tic for the 7 Dor variables, and we propose that HD 49434 
is indeed a variable of this type. A longer time series is 
needed to be able to detect any periodic pulsation signal. 

The star has been chosen as a main target for COROT 
which means that the star will be observed for around 150 
days. We will then be able to probe the interior of this star 
in detail by using seismic techniques. To do this one needs 
to constrain the fundamental parameters of the star which 
indeed has been aim of the work presented here. 
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